Viral Sequences Research Articles

Viral genome sequencing methods: benefits and pitfalls of current approaches.

Whole genome sequencing of viruses provides high-resolution molecular insights, enhancing our understanding of viral genome function and phylogeny. Beyond fundamental research, viral sequencing is increasingly vital for pathogen surveillance, epidemiology, and clinical applications. As sequencing methods rapidly evolve, the diversity of viral genomics applications and catalogued genomes continues to expand. Advances in long-read, single molecule, real-time sequencing methodologies present opportunities to sequence contiguous, haplotype resolved viral genomes in a range of research and applied settings. Here we present an overview of nucleic acid sequencing methods and their applications in studying viral genomes. We emphasise the advantages of different viral sequencing approaches, with a particular focus on the benefits of third-generation sequencing technologies in elucidating viral evolution, transmission networks, and pathogenesis.

Evidence of reassortment of avian influenza A (H2) viruses in Brazilian shorebirds.

Influenza A viruses of the H2 subtype represent a zoonotic and pandemic threat to humans due to a lack of widespread specific immunity. Although A(H2) viruses that circulate in wild bird reservoirs are distinct from the 1957 pandemic A(H2N2) viruses, there is concern that they could impact animal and public health. There is limited information on AIVs in Latin America, and next to nothing about H2 subtypes in Brazil. In the present study, we report the occurrence and genomic sequences of two influenza A viruses isolated from wild-caught white-rumped sandpipers (Calidris fuscicollis). One virus, identified as A(H2N1), was isolated from a bird captured in Restinga de Jurubatiba National Park (PNRJ, Rio de Janeiro), while the other, identified as A(H2N2), was isolated from a bird captured in Lagoa do Peixe National Park (PNLP, Rio Grande do Sul). DNA sequencing and phylogenetic analysis of the obtained sequences revealed that each virus belonged to distinct subtypes. Furthermore, the phylogenetic analysis indicated that the genomic sequence of the A(H2N1) virus isolated from PNRJ was most closely related to other A(H2N1) viruses isolated from North American birds. On the other hand, the A(H2N2) virus genome recovered from the PNLP-captured bird exhibited a more diverse origin, with some sequences closely related to viruses from Iceland and North America, and others showing similarity to virus sequences recovered from birds in South America. Viral genes of diverse origins were identified in one of the viruses, indicating local reassortment. This suggests that the extreme South of Brazil may serve as an environment conducive to reassortment between avian influenza virus lineages from North and South America, potentially contributing to an increase in overall viral diversity.

First report of pepper chlorosis-associated virus infecting tobacco (Nicotiana tabacum) plants in Sichuan Province in China.

Tobacco (Nicotiana tabacum) is an economically important crop in China, and more than 30 viruses have been reported to infect tobacco (Yin et al. 2022). In July 2022, we observed interveinal necrosis on tobacco leaves in fields in Sichuan Province (N 27.9172, E 105.6662) (Fig. 1). Total RNA was isolated from multiple leaves of one plant using an RNAprep Pure Polysaccharide Polyphenol Plant Total RNA Extraction Kit (TIANGEN, Beijing, China). Total RNAs were pooled, and a TruSeq Stranded Total RNA with RiboZero Gold Kit (Illumina, San Diego, CA, USA) was used to eliminate ribosomal RNA. An RNA-Seq library was constructed using VAHTS Universal V6 RNA-seq Library Prep (Nanjing Vazyme, China). High-throughput sequencing was performed on the Illumina DNBseq platform (BGI-ShenZhen, China), which yielded 20,102,087 reads with an average length of 150 nt (total size >6 Gb). Unaligned reads were assembled de novo using SPAdes (Bankevich et al. 2012). Contigs with length ≥200 nt were subjected to local BLASTn and BLASTx analyses against the GenBank nt and nr databases, respectively (Wang et al. 2022). A total of 23 contigs were identified through BLASTx (e-value cut-off = 10 -3), ranging from 631 to 1555 bp long, with 82% to 96% coverage to partial genomic sequences of pepper chlorosis-associated virus (PepCaV-Higashitsuno_2021; Accessions: LC719619 to LC719621) and one contig (6459 bp) with 99% similarity to tobacco mosaic virus (Accession: OP525281) isolate DSMZ PV-0109 from Germany. The complete genome sequence of PepCaV was obtained using primers based on the assembled contigs. The 5'- and 3'-terminal regions of the RNA genome were obtained by 5'- and 3'-rapid amplification of cDNA ends. These amplicons were cloned using the pEASY-Blunt Zero Cloning Kit (TRANSGEN, Nangjing, China) and sequenced by Sanger sequencing. Complete genome sequences of tripartite PepCaV from tobacco samples were 7697, 1808, and 1557 nucleotides long (Accession: OR451987 to OR451989) and showed genome organization typical of the genus Ophiovirus in the family Aspiviridae. The complete sequences of RNA1, RNA2 and RNA3 genome segments shared 92.36%, 90.43%, and 95.24%, nucleotide sequence identities, respectively, with the isolate PepCaV-Higashitsuno_2021 pepper isolate (Accession: LC719619 to LC719621) (Shimomoto et al. 2023), but PepCaV has not been reported to infect N. tabacum. In June 2023, 10 plants collected from each place of Macheng (N 27.9094, E 105.6740), Xiangyang (N 28.0936, E 105.6249) and Moni (N 27.8899, E 105.5936) showing interveinal necrosis symptoms were tested using RT-PCR using PepCaV-MP610-F (5'-TGTTCTCTGCTATGCGGTTG -3') and PepCaV-MP610-R (5'-AGCAATCTCGCACCTGAAGT-3') to product 610bp amplicon. Twenty-five tobacco plants were positive for PepCaV. Single sequence from each location was submitted to GenBank (Accession: PP728631 to PP728633). Sap extracts from the original field leaf samples collected from Sichuan Province were used to mechanically inoculate tobacco plants (10 plants) at the four-leaf stage. After 7 days, leaf samples were tested using RT-PCR assay specific to PepCaV and TMV while samples were positive only for TMV but failed to transmit PepCaV by mechanical inoculation. According to previous literature, ophioviruses may be transmitted though soilborne fungus (Jeong et al. 2014). Further research is needed to understand the transmission, epidemiology, and pathological properties of the PepCaV. To our knowledge, this is the first report documenting natural PepCaV infection of tobacco plants in China, providing a scientific basis for PepCaV infection control in tobacco plantations.

Genome sequence of a vesicular stomatitis Indiana virus isolate collected in 1988 from a naturally infected bovine in Mexico.

We report the near-full genome sequence of a vesicular stomatitis Indiana virus (VSIV) originally collected from a naturally infected bovine in south-central Mexico. This sequence represents a coding-complete genome sequence of a VSIV from Mexico, a country where vesicular stomatitis is endemic.

Characterization of nsp1 Binding to the Viral RNA Leader Sequence of Severe Acute Respiratory Syndrome Coronavirus.

Nonstructural protein 1 (nsp1) of the severe acute respiratory syndrome coronavirus (SCOV1 and SCOV2) acts as a host shutoff protein by blocking the translation of host mRNAs and triggering their decay. Surprisingly, viral RNA, which resembles host mRNAs containing a 5'-cap and a 3'-poly(A) tail, escapes significant translation inhibition and RNA decay, aiding viral propagation. Current literature proposes that, in SCOV2, nsp1 binds the viral RNA leader sequence, and the interaction may serve to distinguish viral RNA from host mRNA. However, a direct binding between SCOV1 nsp1 and the corresponding RNA leader sequence has not been established yet. Here, we show that SCOV1 nsp1 binds to the SCOV1 RNA leader sequence but forms multiple complexes at a high concentration of nsp1. These complexes are marginally different from complexes formed with SCOV2 nsp1. Finally, mutations of the RNA stem-loop did not completely abolish RNA binding by nsp1, suggesting that an RNA secondary structure is more important for binding than the sequence itself. Understanding the nature of binding of nsp1 to viral RNA will allow us to understand how this viral protein selectively suppresses host gene expression.

A highly efficient Agrobacterium-mediated infectious system for Sweet potato leaf curl virus and a deltasatellite-based VIGS vector

Plant virus disease is one of the major threats to the yield and quality of sweetpotato (Ipomoea batatas [L.] Lam.), the seventh most important food crop in the world in terms of production. Sweepoviruses are a group of monopartite geminiviruses (genus Begomovirus, family Geminiviridae) that infect Ipomoea plants and cause significant damage. However, investigations on the mechanisms of sweepovirus infection are scarce due to the low efficiency of the current infectious system. In this work, we used the genomic sequence of Sweet potato leaf curl virus (SPLCV) to construct an effective infectious clone, SPLCV-1.01. We optimized the inoculation procedure in Nicotiana benthamiana by improving the Agrobacterium virulence using an induction medium. Finally, we established a highly efficient agroinoculation protocol for SPLCV on sweetpotato. Furthermore, a simple and reliable “Agro-soaking” method was developed for SPLCV on sweetpotato and Ipomoea aquatica. The “Agro-soaking” method also worked well with Tomato yellow leaf curl China virus (TYLCCNV) and Beet severe curly top virus (BSCTV) on tomato. Based on the “Agro-soaking” method and a deltasatellite, the SBG51 VIGS vector worked with the SPLCV-1.01 infectious clone and silenced two endogenous genes IbCHL (magnesium-chelatase subunit) and IbPDS (phytoene desaturase) in two weeks. Moreover, the VIGS system identified IbADK1 (adenosine kinase 1) as a positive regulator to dampen the SPLCV accumulation in sweetpotato. Our work provides a practical method for studying the interaction between sweetpotato and SPLCV and the functional studies of genes in sweetpotato.

Multimodal data fusion using sparse canonical correlation analysis and cooperative learning: a COVID-19 cohort study

Through technological innovations, patient cohorts can be examined from multiple views with high-dimensional, multiscale biomedical data to classify clinical phenotypes and predict outcomes. Here, we aim to present our approach for analyzing multimodal data using unsupervised and supervised sparse linear methods in a COVID-19 patient cohort. This prospective cohort study of 149 adult patients was conducted in a tertiary care academic center. First, we used sparse canonical correlation analysis (CCA) to identify and quantify relationships across different data modalities, including viral genome sequencing, imaging, clinical data, and laboratory results. Then, we used cooperative learning to predict the clinical outcome of COVID-19 patients: Intensive care unit admission. We show that serum biomarkers representing severe disease and acute phase response correlate with original and wavelet radiomics features in the LLL frequency channel (cor(Xu1, Zv1) = 0.596, p value < 0.001). Among radiomics features, histogram-based first-order features reporting the skewness, kurtosis, and uniformity have the lowest negative, whereas entropy-related features have the highest positive coefficients. Moreover, unsupervised analysis of clinical data and laboratory results gives insights into distinct clinical phenotypes. Leveraging the availability of global viral genome databases, we demonstrate that the Word2Vec natural language processing model can be used for viral genome encoding. It not only separates major SARS-CoV-2 variants but also allows the preservation of phylogenetic relationships among them. Our quadruple model using Word2Vec encoding achieves better prediction results in the supervised task. The model yields area under the curve (AUC) and accuracy values of 0.87 and 0.77, respectively. Our study illustrates that sparse CCA analysis and cooperative learning are powerful techniques for handling high-dimensional, multimodal data to investigate multivariate associations in unsupervised and supervised tasks.

Random-effects substitution models for phylogenetics via scalable gradient approximations.

Phylogenetic and discrete-trait evolutionary inference depend heavily on an appropriate characterization of the underlying character substitution process. In this paper, we present random-effects substitution models that extend common continuous-time Markov chain models into a richer class of processes capable of capturing a wider variety of substitution dynamics. As these random-effects substitution models often require many more parameters than their usual counterparts, inference can be both statistically and computationally challenging. Thus, we also propose an efficient approach to compute an approximation to the gradient of the data likelihood with respect to all unknown substitution model parameters. We demonstrate that this approximate gradient enables scaling of sampling-based inference, namely Bayesian inference via Hamiltonian Monte Carlo, under random-effects substitution models across large trees and state-spaces. Applied to a dataset of 583 SARS-CoV-2 sequences, an HKY model with random-effects shows strong signals of nonreversibility in the substitution process, and posterior predictive model checks clearly show that it is a more adequate model than a reversible model. When analyzing the pattern of phylogeographic spread of 1441 influenza A virus (H3N2) sequences between 14 regions, a random-effects phylogeographic substitution model infers that air travel volume adequately predicts almost all dispersal rates. A random-effects state-dependent substitution model reveals no evidence for an effect of arboreality on the swimming mode in the tree frog subfamily Hylinae. Simulations reveal that random-effects substitution models can accommodate both negligible and radical departures from the underlying base substitution model. We show that our gradient-based inference approach is over an order of magnitude more time efficient than conventional approaches.

Neoantigen Identification and Dendritic Cell-Based Vaccines for Lung Cancer Immunotherapy

Immunotherapies can treat many cancers, including difficult-to-treat cases such as lung cancer. Due to its tolerability, long-lasting therapeutic responses, and efficacy in a wide spectrum of patients, immunotherapy can also help to treat lung cancer, which has few treatment choices. Tumor-specific antigens (TSAs) for cancer vaccinations and T-cell therapies are difficult to discover. Neoantigens (NeoAgs) from genetic mutations, irregular RNA splicing, protein changes, or viral genetic sequences in tumor cells provide a solution. NeoAgs, unlike TSAs, are non-self and can cause an immunological response. Next-generation sequencing (NGS) and bioinformatics can swiftly detect and forecast tumor-specific NeoAgs. Highly immunogenic NeoAgs provide personalized or generalized cancer immunotherapies. Dendritic cells (DCs), which originate and regulate T-cell responses, are widely studied potential immunotherapeutic therapies for lung cancer and other cancers. DC vaccines are stable, reliable, and safe in clinical trials. The purpose of this article is to evaluate the current status, limitations, and prospective clinical applications of DC vaccines, as well as the identification and selection of major histocompatibility complex (MHC) class I and II genes for NeoAgs. Our goal is to explain DC biology and activate DC manipulation to help researchers create extremely potent cancer vaccines for patients.

Ocular manifestations of Monkeypox virus (MPXV) infection with viral persistence in ocular samples – a case series.

ObjectiveWe describe the clinical presentation and ocular viral dynamics in patients with Monkeypox virus-related ophthalmic disease (MPXROD). MethodsIn this case series, we investigated five consecutive patients with confirmed mpox, diagnosed through a positive Monkeypox virus (MPXV) Polymerase Chain Reaction (PCR) test and presenting with ocular symptoms. They were referred from the Reference Center for Sexually Transmitted Infections in São Paulo (CRT) to the Uveitis Sector at the Federal University of São Paulo, between August and December 2022. We performed PCR testing on ocular samples and culture supernatants for MPXV in all patients. Viral sequencing was conducted in one of the cases. ResultsReplicating MPXV was identified in at least one ocular sample of all patients, between day 31 and day 145 after the onset of skin lesions. All patients presented with keratitis, 3 with uveitis (60%) and two exhibited hypopyon (40%). The onset of ocular symptoms occurred at a mean of 21,2 days after the appearance of the first skin lesion and persisted, on average, for 61,6 days, with a worsening trend observed until the initiation of tecovirimat treatment. Tecovirimat treatment was administered to all patients, with initiation occurring between 31 and 145 days after the onset of skin lesions. MPXV genome sequencing of an isolate from one patient classified it as belonging to lineage B1 in clade IIb. ConclusionThis study reveals a late onset and persistence of sight threatening ocular disease, along with potential viral infectivity even after systemic resolution in mpox cases. These findings highlight the risk of ongoing transmission from individuals with prolonged ocular manifestations, particularly through ocular discharge.

Evaluation of the Impact of Concentration and Extraction Methods on the Targeted Sequencing of Human Viruses from Wastewater.

Sequencing human viruses in wastewater is challenging due to their low abundance compared to the total microbial background. This study compared the impact of four virus concentration/extraction methods (Innovaprep, Nanotrap, Promega, and Solids extraction) on probe-capture enrichment for human viruses followed by sequencing. Different concentration/extraction methods yielded distinct virus profiles. Innovaprep ultrafiltration (following solids removal) had the highest sequencing sensitivity and richness, resulting in the successful assembly of several near-complete human virus genomes. However, it was less sensitive in detecting SARS-CoV-2 by digital polymerase chain reaction (dPCR) compared to Promega and Nanotrap. Across all preparation methods, astroviruses and polyomaviruses were the most highly abundant human viruses, and SARS-CoV-2 was rare. These findings suggest that sequencing success can be increased using methods that reduce nontarget nucleic acids in the extract, though the absolute concentration of total extracted nucleic acid, as indicated by Qubit, and targeted viruses, as indicated by dPCR, may not be directly related to targeted sequencing performance. Further, using broadly targeted sequencing panels may capture viral diversity but risks losing signals for specific low-abundance viruses. Overall, this study highlights the importance of aligning wet lab and bioinformatic methods with specific goals when employing probe-capture enrichment for human virus sequencing from wastewater.

Detection, sequencing, and tissue distribution of piscine orthoreovirus 2-like virus in diseased coho salmon in Alaska.

We performed a diagnostic disease investigation on a cohort of coho salmon (Oncorhynchus kisutch) fingerlings in Alaska exhibiting anorexia, gaping mouths, anemia, and increased mortality. Histologic examination revealed mild-to-severe myocardial degeneration and lymphohistiocytic and neutrophilic myocarditis, moderate splenic histiocytosis, and mild renal histiocytosis. Piscine orthoreoviruses 1 and 3 were not detected by molecular methods, and no other viruses could be cultured on 3 common diagnostic fish cell lines. De novo assembly produced a viral genome of 10 linear segments with >80% homology to piscine orthoreovirus 2 (PRV2) encoding all 11 PRV2 proteins. An in situ hybridization probe using RNAscope was developed against 697 viral nucleotides identified by sequencing, which revealed viral genome in heart, spleen, gill, kidney, liver, blood, and the lamina propria of the intestines. Our findings are supportive of a novel piscine orthoreovirus most closely related to PRV2 associated with morbidity and mortality of coho salmon in the northeastern Pacific.

M-F noncoding region sequences of H1 genotype measles virus provide higher resolution for virus transmission tracing

As countries and regions move toward measles elimination, extended sequence window including noncoding region located between the matrix and fusion protein genes (M-F NCR) was considered to be used in molecular surveillance. The molecular resolution of M-F NCR was evaluated with 192 genotype H1 strains circulating during 2011-2018 in China. Phylogenetic analyses of the N450 and M-F NCR targets indicated that both two targets could confirm epi-linked outbreak, while M-F NCR target could further improve resolution of the molecular characterization: (1) it could differentiate the strains with identical N450 circulated in one county within one month of disease onset; (2) different transmission chains could be distinguished for strains with identical N450; (3) better spatial-temporal consistency with topology could be provided among sporadic cases with inconsistent N450. Accordingly, M-F NCR could be used to complement the information from N450 to address the specific questions in tracking the virus transmission chains.

Effect of human H3N2 influenza virus reassortment on influenza incidence and severity during the 2017–18 influenza season in the USA: a retrospective observational genomic analysis

During the 2017-18 influenza season in the USA, there was a high incidence of influenza illness and mortality. However, no apparent antigenic change was identified in the dominant H3N2 viruses, and the severity of the season could not be solely attributed to a vaccine mismatch. We aimed to investigate whether the altered virus properties resulting from gene reassortment were underlying causes of the increased case number and disease severity associated with the 2017-18 influenza season. Samples included were collected from patients with influenza who were prospectively recruited during the 2016-17 and 2017-18 influenza seasons at the Johns Hopkins Hospital Emergency Departments in Baltimore, MD, USA, as well as from archived samples from Johns Hopkins Health System sites. Among 647 recruited patients with influenza A virus infection, 411 patients with whole-genome sequences were available in the Johns Hopkins Center of Excellence for Influenza Research and Surveillance network during the 2016-17 and 2017-18 seasons. Phylogenetic trees were constructed based on viral whole-genome sequences. Representative viral isolates of the two seasons were characterised in immortalised cell lines and human nasal epithelial cell cultures, and patients' demographic data and clinical outcomes were analysed. Unique H3N2 reassortment events were observed, resulting in two predominant strains in the 2017-18 season: HA clade 3C.2a2 and clade 3C.3a, which had novel gene segment constellations containing gene segments from HA clade 3C.2a1 viruses. The reassortant re3C.2a2 viruses replicated with faster kinetics and to a higher peak titre compared with the parental 3C.2a2 and 3C.2a1 viruses (48 h vs 72 h). Furthermore, patients infected with reassortant 3C.2a2 viruses had higher Influenza Severity Scores than patients infected with the parental 3C.2a2 viruses (median 3·00 [IQR 1·00-4·00] vs 1·50 [1·00-2·00]; p=0·018). Our findings suggest that the increased severity of the 2017-18 influenza season was due in part to twointrasubtypes, cocirculating H3N2 reassortant viruses with fitness advantages over the parental viruses. This information could help inform future vaccine development and public health policies. The Center of Excellence for Influenza Research and Response in the US, National Science and Technology Council, and Chang Gung Memorial Hospital in Taiwan.

Identification and epidemiological study of an uncultured flavivirus from ticks using viral metagenomics and pseudoinfectious viral particles

During their blood-feeding process, ticks are known to transmit various viruses to vertebrates, including humans. Recent viral metagenomic analyses using next-generation sequencing (NGS) have revealed that blood-feeding arthropods like ticks harbor a large diversity of viruses. However, many of these viruses have not been isolated or cultured, and their basic characteristics remain unknown. This study aimed to present the identification of a difficult-to-culture virus in ticks using NGS and to understand its epidemic dynamics using molecular biology techniques. During routine tick-borne virus surveillance in Japan, an unknown flaviviral sequence was detected via virome analysis of host-questing ticks. Similar viral sequences have been detected in the sera of sika deer and wild boars in Japan, and this virus was tentatively named the Saruyama virus (SAYAV). Because SAYAV did not propagate in any cultured cells tested, single-round infectious virus particles (SRIP) were generated based on its structural protein gene sequence utilizing a yellow fever virus-based replicon system to understand its nationwide endemic status. Seroepidemiological studies using SRIP as antigens have demonstrated the presence of neutralizing antibodies against SAYAV in sika deer and wild boar captured at several locations in Japan, suggesting that SAYAV is endemic throughout Japan. Phylogenetic analyses have revealed that SAYAV forms a sister clade with the Orthoflavivirus genus, which includes important mosquito- and tick-borne pathogenic viruses. This shows that SAYAV evolved into a lineage independent of the known orthoflaviviruses. This study demonstrates a unique approach for understanding the epidemiology of uncultured viruses by combining viral metagenomics and pseudoinfectious viral particles.

Performance of amplicon and capture based next-generation sequencing approaches for the epidemiological surveillance of Omicron SARS-CoV-2 and other variants of concern.

To control the SARS-CoV-2 pandemic, healthcare systems have focused on ramping up their capacity for epidemiological surveillance through viral whole genome sequencing. In this paper, we tested the performance of two protocols of SARS-CoV-2 nucleic acid enrichment, an amplicon enrichment using different versions of the ARTIC primer panel and a hybrid-capture method using KAPA RNA Hypercap. We focused on the challenge of the Omicron variant sequencing, the advantages of automated library preparation and the influence of the bioinformatic analysis in the final consensus sequence. All 94 samples were sequenced using Illumina iSeq 100 and analysed with two bioinformatic pipelines: a custom-made pipeline and an Illumina-owned pipeline. We were unsuccessful in sequencing six samples using the capture enrichment due to low reads. On the other hand, amplicon dropout and mispriming caused the loss of mutation G21987A and the erroneous addition of mutation T15521A respectively using amplicon enrichment. Overall, we found high sequence agreement regardless of method of enrichment, bioinformatic pipeline or the use of automation for library preparation in eight different SARS-CoV-2 variants. Automation and the use of a simple app for bioinformatic analysis can simplify the genotyping process, making it available for more diagnostic facilities and increasing global vigilance.

First report of cucumber green mottle mosaic virus infecting Cynanchum rostellatum in China.

Cucumber green mottle mosaic virus (CGMMV) was first discovered on cucumber in the United Kingdom in 1935 (Ainsworth, 1935), and has spread worldwide except to Antarctica (Jones, 2021). Given its extensive damage, it is considered an important pathogen on global cucurbit plants and fruit crops. In China, CGMMV was first reported on pumpkin in Guangxi Province in 2003 (Qin et al., 2005), and occurred on 34 plants species across 23 provinces (Liu et al., 2016). Cynanchum rostellatum is a member of the family Apocynaceae. In July 2021, leaves of C. rostellatum exhibiting virus-like symptoms (yellowing, severe crinkling, deformation) were observed and collected in Liaoning Province, China. Aphids were also observed on the leaves and stems (Fig. S1) of the plants and were collected. Total RNA was extracted from diseased leaves following the CTAB method, followed by the depletion of ribosomal RNAs (rRNA) with TIANSeq rRNA Depletion Kit (Tiangen, China). The RNAs were, then processed into a DNBSEQ LncRNA-Seq library, and sequenced on the MGISEQ-2000 platform at BGI Genomics (Wuhan, China). A total of 106.98 M clean reads were obtained after data filtering using SOAPnuke software (BGI, China). The clean reads were assembled into contigs using CLC Genomics Workbench 11 (Qiagen, USA) and Trinity v2.0.6 (Haas et al., 2013). A contig (4,760 reads, average coverage:73.76) of 6,391 nucleotides was found to share the highest sequence identity (99.83%) with CGMMV isolate GDLZ (MK933286), irrespective of other virus-like contigs related to Polerovirus and Totivirus. Based on the genome of GDLZ isolate, seven specific primers (Table S1) were designed to amplify the full viral genomic sequences using a PrimeScriptTM One-Step RT-PCR Kit. Seven expected amplicons were obtained, cloned, and sequenced. The complete genome was determined to be 6,423 nucleotides (GenBank accession number OR854819) in length and designated as LNMJ isolate. LNMJ shared 96.8%-99.7% nucleotide sequence identities with CGMMV isolates from China. Phylogenetic analysis based on the complete genome sequences showed that LNMJ clustered together with CGMMV isolates hn (GenBank accession number KC851866), GDLZ (GenBank accession number MK933286), and JD8 (GenBank accession number KM873784) from China. The specific primers LM-TJ-3F/3R were designed to determine the virus-symptom association for LNMJ, and all twelve symptomatic C. rostellatum plants collected from fields tested positive for LNMJ. Two out of six randomly selected aphids from the diseased plants also tested positive. To further prove its infectivity, LNMJ was inoculated mechanically onto ten healthy Nicotiana benthamiana plants, and the results indicated a high infection rate of 80% (8/10), at 30 days post-inoculation despite no distinct symptoms observed. To our knowledge, this is the first report of the natural infection of C. rostellatum plants with CGMMV. C. rostellatum is a widespread herb in China (Wei et al., 2019) and more surveys are needed to determine the distribution of CGMMV. The habitats of C. rostellatum span diverse agroecological zones, and thus our study underscores the potential spillover of CGMMV to neighboring crops as a significant risk.

A hepatitis B virus (HBV) sequence variation graph improves alignment and sample-specific consensus sequence construction.

Nearly 300 million individuals live with chronic hepatitis B virus (HBV) infection (CHB), for which no curative therapy is available. As viral diversity is associated with pathogenesis and immunological control of infection, improved methods to characterize this diversity could aid drug development efforts. Conventionally, viral sequencing data are mapped/aligned to a reference genome, and only the aligned sequences are retained for analysis. Thus, reference selection is critical, yet selecting the most representative reference a priori remains difficult. We investigate an alternative pangenome approach which can combine multiple reference sequences into a graph which can be used during alignment. Using simulated short-read sequencing data generated from publicly available HBV genomes and real sequencing data from an individual living with CHB, we demonstrate alignment to a phylogenetically representative 'genome graph' can improve alignment, avoid issues of reference ambiguity, and facilitate the construction of sample-specific consensus sequences more genetically similar to the individual's infection. Graph-based methods can, therefore, improve efforts to characterize the genetics of viral pathogens, including HBV, and have broader implications in host-pathogen research.

Comparing the molecular evolution and recombination patterns of predominant PRRSV-2 lineages co-circulating in China.

Porcine reproductive and respiratory syndrome virus (PRRSV) poses widespread epidemics in swine herds, yet the drivers underlying lineage replacements/fitness dynamics remain unclear. To delineate the evolutionary trajectories of PRRSV-2 lineages prevalent in China, we performed a comprehensive longitudinal phylodynamic analysis of 822 viral sequences spanning 1991-2022. The objectives encompassed evaluating lineage dynamics, genetic diversity, recombination patterns and glycosylation profiles. A significant shift in the dominance of PRRSV-2 sub-lineages has been observed over the past 3 decades, transitioning from sub-lineage 8.7 to sub-lineage 1.8, followed by extensive diversification. The analysis revealed discordant recombination patterns between the two dominant viral sub-lineages 1.8 and 8.7, underscoring that modular genetic exchanges contribute significantly to their evolutionary shaping. Additionally, a strong association was found between recombination breakpoint locations and transcriptional regulatory sequences (TRSs). Glycosylation patterns also demonstrated considerable variability across sub-lineages and temporally, providing evidence for immune-driven viral evolution. Furthermore, we quantified different evolutionary rates across sub-lineages, with sub-lineage 1.8 uniquely displaying the highest nucleotide substitution rates. Taken together, these findings provide refined insight into the evolutionary mechanisms underpinning cyclic shifts in dominance among regionally circulating PRRSV sub-lineages.

COVID-19 cluster surveillance using exposure data collected from routine contact tracing: The genomic validation of a novel informatics-based approach to outbreak detection in England.

Contact tracing was used globally to prevent onwards transmission of COVID-19. Tracing contacts alone is unlikely to be sufficient in controlling community transmission, due to the pre-symptomatic, overdispersed and airborne nature of COVID-19 transmission. We describe and demonstrate the validity of a national enhanced contact tracing programme for COVID-19 cluster surveillance in England. Data on cases occurring between October 2020 and September 2021 were extracted from the national contact tracing system. Exposure clusters were identified algorithmically by matching ≥2 cases attending the same event, identified by matching postcode and event category within a 7-day rolling window. Genetic validity was defined as exposure clusters with ≥2 cases from different households with identical viral sequences. Exposure clusters were fuzzy matched to the national incident management system (HPZone) by postcode and setting description. Multivariable logistic regression modelling was used to determine cluster characteristics associated with genetic validity. Over a quarter of a million (269,470) exposure clusters were identified. Of the eligible clusters, 25% (3,306/13,008) were genetically valid. 81% (2684/3306) of these were not recorded on HPZone and were identified on average of one day earlier than incidents recorded on HPZone. Multivariable analysis demonstrated that exposure clusters occurring in workplaces (aOR = 5·10, 95% CI 4·23-6·17) and education (aOR = 3·72, 95% CI 3·08-4·49) settings were those most strongly associated with genetic validity. Cluster surveillance using enhanced contact tracing in England was a timely, comprehensive and systematic approach to the detection of transmission events occurring in community settings. Cluster surveillance can provide intelligence to stakeholders to support the assessment and management of clusters of COVID-19 at a local, regional, and national level. Future systems should include predictive modelling and network analysis to support risk assessment of exposure clusters to improve the effectiveness of enhanced contract tracing for outbreak detection.